Topological Bose-Mott Insulators in a One-Dimensional Optical Superlattice 
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We study topological properties of the Bose-Hubbard model with repulsive interactions in a one- 
dimensional optical superlattice. We find that the Mott insulator states of the single-component 
(two-component) Bose-Hubbard model under fractional fillings are topological insulators character- 
ized by a nonzero charge (or spin) Chern number with nontrivial edge states. For ultracold atomic 
experiments, we show that the topological Chern number can be detected through measuring the 
density profiles of the bosonic atoms in a harmonic trap. 
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Introduction - Ultracold atoms in optical lattices can 
be used to simulate strongly correlated many-body mod- 
els that are central to the understanding of condensed 
matter physics. This simulation has attracted a lot 
of attention as the optical lattice experiments offer un- 
paralleled controllability and new tools to study many- 
body physics [3-0|. As a remarkable example, the Bose- 
Hubbard (BH) model has been experimentally realized 
with ultracold atoms and a quantum phase transition 
from a superfluid to a Mott insulator described by this 
model has been observed [3] . On the other hand, topolog- 
ical matters, such as quantum Hall systems and topolog- 
ical insulators, are of fundamental importance in physics 
[B| . Recently, studying topological phases with ultracold 
atoms has raised great interest 0, H-12|. In general, it 
requires complicated control of experimental systems to 
realize topological phases with ultracold atoms. An in- 
teresting question is whether one can observe topological 
phases in a simple BH type of model, which can be read- 
ily implemented by many experimental groups. Topo- 
logical properties of bosonic systems, however, have not 
been well-studied in literature, partly for the reason that 
the topological invariants are usually defined as an in- 
tegration over all the occupied states in the momentum 
space 13|, llj] . This definition does not apply directly to 
the bosonic system as many bosons can occupy the same 
momentum state. 

In this Letter, in contrast to the conventional wisdom, 
we show that the BH model in a one-dimensional (ID) 
optical superlattice displays nontrivial topological prop- 
erties. We demonstrate that the Mott insulators of the 
single-component (two-component) BH model at frac- 
tional fillings belong to topological matter with its phase 
characterized by a nonzero integer charge (or spin) Chern 
number. For Mott insulators, the bulk system has a gap 
in the excitation spectrum induced by the interaction. 



For a topologically nontrivial Mott insulator state charac- 
terized by a nonzero Chern number, we further show that 
there are protected edge states inside the bulk gap under 
the open boundary condition. The Mott insulators at in- 
teger fillings for this system remain topologically trivial 
with a zero Chern number and no edge states. The topo- 
logical properties discussed here are reminiscent of those 
in topological Mott insulators theoretically predicted in 
Ref. 15|] for the Fermi-Hubbard model in a honeycomb 
lattice with frustrated next-neighbor interactions. Re- 
markably, we here show that topological Mott insulators 
can appear in a simple ID BH model in an optical su- 
perlattice, which, besides being conceptually interesting, 
makes the experimental realization of topological matter 
much easier in the ultracold atomic system. We pro- 
pose a scheme to detect the topological Chern number 
by observation of the plateaus of the density profile with 
ultracold atoms in a weak global harmonic trap as it is 
the case for experiments. 

Single-component BH model in a superlattice.- We 
consider a single-component bosonic gas loaded into a 
ID optical lattice, which is described by the BH model 



H 



+ J2[Un j (n j -l)/2 + V j n j ] 
3 



(1) 



where V 3 ■ — V cos(2irctj + 8) denotes a periodic superlat- 
tice potential [l6| . bj and &j correspond to the bosonic 
annihilation and creation operators of atoms on the jth 
lattice site, nj — bjbj is the number operator, and J and 
U represent the hopping rate and the on site interaction 
strength, respectively. We consider in this Letter a com- 
mensurate superlattice potential Vj with a = p/q (p,q 
are integers) being a rational number and 8 an arbitrary 
tunable phase, which has been experimentally realized 
16]. We take J as the energy unit by setting J = 1. 
The ground-state phase diagram of the Hamiltonian 



2 



[Eq.QJ] is well-studied di, Il7j ]. For a sufficiently large 
[/, the system is in a gapped Mott insulator phase at 
commensurate fillings with v = N^/N = ma, where m is 
an integer, Nf, is the atom number, and N is the number 
of lattice sites. Away from the commensurate fillings or 
for a small U, the system is in a superfluid state [17|. In 
this Letter, we focus on study of the topological proper- 
ties of the Mott insulator phase. 

The energy gap and the Chern number of the ground 
state.- The topological property is best characterized by 
the Chern number. To calculate the Chern number, we 
first perform exact diagonalization of the Hamiltonian 
[Eq.JT])] on a chain of N sites with periodic or open 
boundary conditions 



18] 



The ground state is nonde- 
generate and separated from the higher eigenstates by a 
finite gap A at the commensurate fillings. This gap is 
shown in Fig.l as a function of the interaction strength 
U at v — 1/3. The gap increases monotonically with U 
and then saturates at a finite value. For a large U, the 
atoms become hard-core bosons. In this case, each site 
is occupied by no more than one atom. The hard-core 
boson Hubbard model can be mapped to a model of free 
fcrmions. From that mapping, we find that the satura- 
tion value of the gap is 1.08 (in units of J) at a large U 
for an infinite system. The gap should decrease to zero 
as U drops below a critical value U c where the system 
transits to a superfluid phase. In Fig. 1, due to the fi- 
nite size effect, the saturation value of the energy gap is 
above 1.08 for the large U and does not drop exactly to 
zero as U diminishes. However, as the number of lattice 
sites increases, we clearly see the tendency that the gap 
approaches these limiting values at the two ends. 
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FIG. 1: (Color online) The energy gap A and the Chern num- 
ber C 3 [ defined by Eq.(f3])] as functions of the interaction 
strength U. The number of lattice sites N = 9, 12, 15 are 
used for exact diagonalization, and we take TV = 15 for cal- 
culation of C g . Other parameters include V = 1.5, S = 2n/3, 
a — 1/3, and v — 1/3. 



Now, we investigate the topological property of the 
system by calculating the Chern number. For fermions, 
the Chern number is defined as an integration over the 
occupied states in the momentum space [13| . This defini- 
tion can not be extended to the bosonic system as many 
bosons can occupy the same momentum state. Fortu- 



nately, there is another way to calculate the Chern num- 
ber for interacting systems [l9[ : suppose the ground state 
has a gap to the excited state and it depends on the pa- 
rameters 9, S through a generalized periodic boundary 
condition |*(j + N,9,S)) = e id \$(j,6,8)), where j de- 
notes an arbitrary site, 9 is the twist angle, and S is the 
phase in the superlattice potential Vj . Under this bound- 
ary condition, we numerically diagonalize Hamiltonian 
[Eq.JTJ)] and derive the ground state \^(9, 5)), which is a 
non-degenerate state separated from the excited state by 
a nonzero energy gap A when U > U c . For the ground 
state \^f(6,6)) where 9 and 6 vary on a torus, one can 
define the Chern number C g as a topological invariant 
by the following formula [19( 



C -- 
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d9 dS(d e A s -d s A e ), (2) 
Jo 



where the Berry connection = i(^f(6,6)\dfj,\^f(6,6)) 
(fi = 5,9). We numerically calculate the Chern number 
C g using the method for a discrete manifold [2(|. When 
the parameter a — 1/3, we find for this boson system 
that the Chern number C g = 1 (—1) for the filling frac- 
tion v = 1/3 (2/3) and C g = when v = 1. As an 
example, we show the value of C g as a function of U at 
v = 1/3 in Fig. 1, where the manifold of torus is dis- 
cretized by 5 x 5 meshes in the calculation. When the 
system is in the gapped Mott state with U > U c , C g 
is quantized to be exactly at 1, while C' g is unquantized 
when the system enters the gapless superfluid phase. Be- 
cause of quantization of C gi the finite size effect seems 
to have a minimal influence, and we can use exact di- 
agonalization of a small system to get the exact value 
of C g in Fig. 1 for the Mott phase (however, because 
of the finite-size gap, C g is still approximately unity in 
some region of the superfluid phase near the transition 
point). This calculation unambiguously shows that this 
bosonic system is in a topological Mott insulator phase 
with nonzero Chern number at the fractional filling of 
the optical lattice. 

Edge states.- The appearance of edge states at the 
boundary is usually considered to be a hallmark of non- 
trivial topological properties for the bulk system. Under 
the periodic boundary condition, this interacting system 
is gapped at the fractional filling v = 1/3 (or 2/3). How- 
ever, under the open boundary condition, edge states 
confined to the boundary can appear inside the energy 
gap, signaling the nontrivial topological properties of the 
bulk insulator. The quasiparticle energy spectrum AE n 
is determined by the additional energy required to add 
an atom to a system with n atoms, that is, 



A#- p > = E. 



(O.P) 



(3) 



where E^ (E^ ) is the ground-state energy of the sys- 
tem with n atoms under the open (periodic) boundary 
condition [2l|. In Fig. 2(a) , we show the quasi-particle 
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energy spectrum for a system with 96 lattice sites near 
the filling v — 1/3 under both periodic and open bound- 
ary conditions. The calculation is done using the density 
matrix renormalization group method 22J], which pro- 
vides a reliable approach to precisely calculate energies 
for any ID systems. Near the filling v — 1/3, the quasi- 
particle energy spectrum is split into two branches sepa- 
rated by a finite gap. The calculation clearly shows that 
two states appear in the gap of the energy spectrum un- 
der the open boundary condition. In Fig. 2 (b) , we show 
the quasi-particle energy spectrum as a function of phase 
8 under the open boundary condition. Inside the gap be- 
tween the lower and the upper branches of the energy 
spectrum, one can see two edge modes (which are de- 
generate in energy at 5 = 2tt/3) that connect these two 
branches of the bulk spectrum as one varies phase 8. 

To verify that the in-gap modes indeed correspond to 
the edge states, we numerically calculate the excitation 
distribution of these modes and find that they are con- 
fined near the edges of the chain. The distribution of the 
quasi-particle can be defined as 



Arij = 



(4) 



where |^) denotes the ground state wave function of the 
system with n bosonic atoms. The distribution of the in- 
gap quasi-particle modes for N = 96 sites under filling 
v = 1/3 is plotted in Fig. 2 (c). As expected, the in-gap 
states mainly distribute near the two edges, especially for 
a large V. For instance, 99% of the quasi-particle modes 
at V = 10 are localized at the two edge sites. 

Two-component BH model in a superlattice and spin 
Chern number.- If the phase S in Eq.(2) is replaced by 
— 8, we find that the Chern number is —1 (1) for v = 1/3 
(2/3), that is, the sign of the Chern number is flipped. 
This fact implies that we can realize a topological insula- 
tor characterized by a nontrivial spin Chern number with 
a two-component bosonic gas in a ID optical superlattice. 
To this end, we consider a simple case where the inter- 
component atomic collision is turned off, e.g., through 
Feshbach resonance, and the system is described by a 
decoupled two-component BH model with the Hamilto- 
nian 



where the potential V ]a — V^cos(27raj + 8 a ) with Sf = 
—8^ = 8, bj a denotes the bosonic annihilation opera- 
tor with the pseudo-spin a =y, and U a is the intra- 
component interaction rate for spin cr. For this decoupled 
two-component BH model, we have the Chern number 
= — Cj. So, although the total charge Chern num- 
ber Cj + Cjj cancels out to zero, the spin Chern number 
Cg = Cj — Cjj [23j is non- vanishing at fractional fillings. 
The nonzero spin Chern number is usually associated 
with the quantum spin Hall effects in two-dimensional 
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FIG. 2: (Color online) (a) The quasi-particle energy spec- 
trum AE„ [see the definition by Eq. (3)] versus n under the 
periodic (PBC) or open (OBC) boundary condition. The cal- 
culation is done in a 96-site lattice near the filling v = 1/3 
with V = 1.5, U = 10, 5 = 2tt/3, and a = 1/3. (b) The edges 
of the lower (AS30) and the upper (AE33) branches of the 
energy spectrum and the two in-gap modes (Ai^iand A.E32) 
as functions of the phase r5 under the open boundary condi- 
tion, (c) The distribution of the two in-gap modes along the 
chain at V = 1.5, 10. The other parameters for (b) and (c) 
are the same as those for (a). 



systems [23j . For our ID system, spin edge states appear 
when the spin Chern number is nonzero. For an example 
with a = 1/3, we have Cg — 2 at the fractional fillings 
v = l/3 (2/3) and C s g = at the integer filling. The edge 
states are similar to those shown in Fig. 2. The spin up 
(down) edge state is confined near the left (right) edge, 
respectively. 

Experimental detection.- We now discuss how to mea- 
sure topological Chern number in a practical experimen- 
tal setting. For atomic experiments, apart from the op- 
tical superlattice potential, the bosons are confined in 
a weak global harmonic trap. For simplicity, we con- 
sider the large- U limit where the system is described by 
hard-core bosons with no more than one atom occupy- 
ing the same lattice site. The total potential, including 
the optical superlattice and the global harmonic trap, is 
described by 



Vj = V cos(27raj + 8) + V H (j - jo) 2 



(6) 



where jo denotes the position of the trap center and 
Vh is the strength of the harmonic trap. We use the 



Jordan- Wigner transformation, b 



3-1 p -inflf„ 
m—1 



and bj = Ilm=i e l7T ^^ m fj, to map the hard-core BH 
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model to non-interacting fermion Hamiltonian Hp = 
■I V, :/;/,., + h.c.) + Ei^/J/i, where /] and ./' 
are the creation and annihilation operators for spinless 
fermions, respectively [tR I24] ]. The particle density of 
hard-core bosons coincides with that of non-interacting 
fermions as we have nj = (bjbj) = (fjfj) = with the 
Jordan- Wigner transformation; however, the momentum 
distribution for bosons is typically very different from 
that for fermions. 

After the hard-core bosons are mapped to fermions, 
there is a simple way to figure out the Chern number. 
The ground state of free fermions is a Slater determi- 
nant, i.e., a product of single particle states |\&^) = 



1/3,1/4, with the results shown in Fig. 3. To reduce 
the oscillations in density profiles induced by modula- 
tion of the potentials, we define the local average den- 
sity fij = J2 m =-M n j+m/ (2M + 1), where 2M + 1 is the 
length to average the density, which corresponds to the 
position resolution in the experimental detection. We 
take M <C N, e.g., M = 4 and N = 300 in Fig. 3, 
as it is typical for experiments. As one can see from 
the density profiles fij in Fig. 3, plateaus appear at 
the rational fillings v = 1/3,2/3,1 for a = 1/3, and 
v = 1/4,3/4,1 for a = 1/4 (the gap at half filling in 
the case of a = p/q with an even q is generally closed 
at an integer S/n [L4j). Using the Streda formula [Eq. 



]lm=i En=i p n m fl\0), with N f = N b the number of ©], we obtain C g = 5n 3 = 1, -1 at the fractional fillings 



fermions and P the matrix of the components of l^g). 
Supposing that the n-th eigenstate of a single particle 
is denoted by \ip n ) = J2j <fij,nfj\0), the eigenvalue equa- 
tion Hf\iP„,) = E n \ip n ) can be written in terms of the 
following Harper equation fllj ] 



J(<l>j+i,n + <f>j,n) + Vcos(2Tra + 5)<f>j tn = E n 4> jt „, (7) 



where 



is the amplitude of the particle wave func- 



tion of the j-th site and E n is the n-th single-particle 
eigen-energy. Compared with the Harper equation in a 
magnetic field, we know that a corresponds to the mag- 
netic flux 14j. Therefore, we can define the local density 
difference as 



drij = 



rijOi) -n j {a 2 ) 
ai - a 2 



(8) 



The Chern number C g can then be obtained through the 
Streda formula C g — Srij under the condition that nj(a^) 
(r) = 1,2) is the local density associated with the plateau 
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FIG. 3: (Color online) The average density profiles n, and 
the density difference Sfij for ai = 1/3 and c<2 = 1/4. The 
values of Sfij represent the corresponding Chern numbers at 
the plateaus with the fillings v — 1/3, 2/3, 1 for a = 1/3 and 
v — 1/4, 3/4, 1 for a = 1/4. The other parameters are iV — 
300, N b = 180, M = 4, V = 10, <5 = tt/2, and Vh = 0.001. 



Following the method outlined in Ref. [T/J, we nu- 



v = a, 1 — a, and C g — at the integer filling v = 1. 
The width of the plateaus is associated with the size of 
the energy gap. To make detection of the Chern number 
easier, we can adjust the frequency of the harmonic trap 
to move the target plateaus to the center of the trap. 
For example, if we choose Vh — 10 -4 J and other param- 
eters as those given in Fig. 3, the plateaus at v = 2/3 for 
a = 1/3 and ^ = 3/4forev = l/4 are moved to the center 
of the trap spanning from the 65th to the 235th lattice 
site. With such a wide plateau, it is straightforward to 
read out the Chern number C g = Sfij = — 1 for this case. 

In summary, we have shown that for bosonic atoms 
in a ID optical super-lattice, the Mott insulator states 
of the corresponding BH model at fractional fillings are 
topologically nontrivial, characterized by nonzero Chern 
number and existence of edge states. We further pre- 
dict that the topological Chern number can be detected 
by measuring the plateaus in the density profile when 
the atoms are trapped in a global harmonic potential. 
The model discussed in this Letter represents one of the 
simplest experimental systems to show intriguing topo- 
logical properties, and the proposed detection method 
allows one to confirm these topological properties with 
the state-of-the-art technology. 
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